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The heat capacity of pure ^He in low density aerogel is measured at 22.5 bar. The superfluid 
response is simultaneously monitored with a torsional oscillator. A slightly rounded heat capacity 
peak, 65 fiK in width, is observed at the ^He-aerogel superfluid transition. Tea- Subtracting the 
bulk '^He contribution, the heat capacity shows a Fermi-liquid form above Tea- The heat capacity 
attributed to superfluid within the aerogel can be flt with a rounded BCS form, and accounts for 
0.30 of the non-bulk fluid in the aerogel, indicating a substantial reduction in the superfluid order 
parameter consistent with earlier superfluid density measurements. 
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The 1995 discovery of the superfluid transition of '^He 
in aerogel glass has aroused wide interest. The 

■^He-aerogel system presents a unique opportunity for 
study of the influence of quenched impurities on the well- 
understood system of pure superfluid '^He. Earlier studies 
employing other porous media such as sintered silver or 
packed powders Q Q were dominated by surface scatter- 
ing and finite size effects ^ . Aerogel consists of a fractal 
structure of relatively uniform silica strands 30 A in di- 
ameter. In the case of 97.6% open aerogel, employed 
in the present experiment, scattering due to the silica 
strands, though strong enough to introduce a significant 
quasiparticle density at T=0, is sufficiently weak so that 
superfluidity is not as strongly suppressed as it would be 
the case of more dense aerogel systems or packed pow- 
ders. While previous measurements exposed the onset of 
superfluidity using torsional pendule m i| NMR g @] 
lH 111 ^^^d acoustic techniques [^| |Q, they could 
not distinguish between a percolation transition, where 
the flow and the specific heat anomaly are distinct, and 
a suppressed superfluid transition p3[ |. 

In the work presented here, we examined the "^He B 
phase superfluid state in aerogel through heat capacity 
measurements at a pressure of 22.5 bar [^,^. Our chief 
result is the observation of a feature in the heat capacity 
coincident with the onset of superflow in the aerogel-'^He 
system. This observation confirms that superflow on- 
set corresponds to a true phase transition rather than to 
a possible dynamic effect associated with a percolation 
transition involving distributed superfluid regions. The 
analysis reveals a normal fermi liquid like contribution to 
the heat capacity down to the lowest temperatures, cor- 
responding to the existence of a non-zero normal fraction 
down to absolute zero. In earlier work [l^ , we reported 
on measurements of the heat capacity employing a drift 
technique. Those measurements gave an indication of an 
anomaly at the aerogel superfluid transition; however, 
a definitive resolution of this peak has required the im- 
provements of the present experiment. 
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FIG. 1. A cross-section of the experimental cell is shown. 

In Fig. 1, we show a cut-away diagram of the latest 
version of our heat capacity apparatus. The cell contain- 
ing the ■^He-aerogel sample also forms the inertial head of 
a torsional oscillator. We are thus able to correlate fea- 
tures in the heat capacity with the dynamic superfluid 
response. The liquid "^He within the apparatus is cooled 
by a silver sinter pad, which is thermally clamped to a 
PrNis nuclear cooling stage. There are two lanthanum 
diluted cerium magnesium nitrate (LCMN) thermome- 
ters; one inside the sample cell and the other located 
near the silver sinter. These are operated in a dc mode 
similar to that of the SQUID based thermometer devel- 
oped by Lipa and Chui IT^ ] . Thermometer calibration is 
provided by comparison to a melting curve thermometer 
mounted on the nuclear cooling stage Jl6[ | . Thermal con- 
tact between the heat capacity cell and the cooling stage 
is established via the '^He contained in the hollow torsion 
rod that connects the heat capacity cell to the stage. The 
time constant associated with this thermal contact path 
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FIG. 2. The thermal response of the cell thermometer to a 
series of four heat pulses (depicted at the bottom of the figure) 
is shown as a function of time. The solid line indicates the 
steady drift of the stage temperature and the dashed curves 
are fits of the exponential decay function (see text) to the 
relaxation portion of the cell's thermal response following a 
heat pulse. 



to restricting these times to the exponential portion 
of the thermal relaxation of the cell. During this period 
the LCMN thermometer gives a reliable value for the 
temperature throughout the cell. The energy flow from 
the cell, due to thermal conduction down the torsion rod 
to the stage, is obtained by numerical integration of the 
quantity K{Tc{t) — Ts{t)) , over the chosen time interval. If 
the time interval contains the heat pulse, we include the 
energy, Qm, contributed by the heater. A possible choice 
for the first time interval would be to take ti = 0, the 
time just before the heat pulse, and t2 = tf, the time at 
the end of the relaxation period just before the next heat 
pulse. The second interval, which will not include a heat 
pulse, might start at time t (chosen in the exponential 
decay region) and ends at tf as for the first interval. 
Then we write two independent equations containing the 
heat capacity and the thermal conductivity. 



C(T,{tf ) ~ r,(0)) = Q,n " / {Tc{t) - T,{t))dt (1) 



is designed to be at least an order of magnitude longer 
than the thermal relaxation time within the heat capac- 
ity cell itself. Heat pulses are applied to the sample by a 
non-inductive resistive heater wound around the aerogel 
sample. Care is taken that the heater is not close to the 
LCMN thermometer. In addition, a second spare heater, 
shown in Fig. 1, was placed in the cell. 

Our measurement technique is illustrated in Fig. 2. 
After a demagnetization to our lowest temperatures, the 
cell warms slowly over several days. During this warm- 
ing interval, a sequence of heat pulses is applied to the 
•^He within the cell and the heat capacity is determined 
as a function of temperature. Representative data are 
shown in Fig. 2, illustrating a series of heat pulses and 
the subsequent thermal relaxation of the sample cell tem- 
perature, Tc(t). In this example, the stage temperature 
Ts{t), is increasing at a steady rate of 15.6 /zK/hr and 
identical heat pulses are repeated at 1800 second inter- 
vals. Following each heat pulse, the cell temperature rises 
rapidly. Initially the temperature distribution within the 
cell is nonuniform. The internal equilibrium time is rel- 
atively short, so that over time the thermal distribution 
within the cell becomes nearly uniform. The cell tem- 
perature exhibits an exponential relaxation of the form, 
Tc{t) = Tqc^*'^'^ +Ts{t) — aT , where a is the warming rate 
of the stage and t is the relaxation time constant. Fits 
of this function to the relaxation portion of the thermal 
pulses are shown in the figure. 

The temperature excursion during each heat pulse is 
relatively small, on the order of 40 to 50 /uK. Therefore, 
in our analysis we take the heat capacity, C, and thermal 
conductivity, k, to be constant over the period of each 
relaxation. We proceed by calculating the net energy 
flow, Q, into or out of the cell during a time interval, ti 



C{T,{tf) - T,{t)) = ^kJ {T,{t) ~ T,{t))dt (2) 

These equations are then solved for C and k with 
A defined as J^{Tc{t) - Ts{t))dt and Af defined as 

J^^{Ut)-Ut))dt. 



c 



(Af - A)Q 



Afinit) - r,(0)) - A{T,itf) ~ T,(0)) 



(T,(t) - T,itf))Q 



Af{TS) - TM) - A{Tc{tf) - T,(0)) 



(3) 



(4) 



One advantage of this approach is that one can vary i, 
and check that the calculated values for C and k are in- 
dependent of t. The thermal conductivity obtained from 
the above analysis is found to be a smoothly varying func- 
tion of temperature except at the bulk superfluid transi- 
tion where the conductivity drops by over a factor of two 
as the temperature rises through the transition tempera- 
ture. A more detailed discussion of the data analysis and 
the cell construction will be found elsewhere |l7j fl^ . 

Throughout the discussion that follows the bulk con- 
tribution is taken to arise from an equivalent amount of 
superfluid in the absence of the aerogel. 

In Fig. 3A, we have plotted the total heat capacity 
determined for temperatures between 0.6 to 3 mK. We 
also show the Q^^ for the torsional oscillator. There are 
two conspicuous features in the heat capacity data. The 
first is the heat capacity anomaly associated with the 
superfluid transition in the '^He-aerogel sample, which 
coincides with the superfiuid aerogel transition, Tca^ as 
marked by the torsional oscillator. The second feature 
is the sharp jump in the heat capacity due to bulk ■^He 
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FIG. 3. Panel A shows the total measured heat capacity for 
the '^He-aerogel sample at a pressure of 22.5 bar, the calcu- 
lated bulk ^He contribution, and the dissipation, Q^^, of the 
torsional oscillator. Panel B shows the heat capacity after the 
bulk contribution is subtracted, along with the BCS heat ca- 
pacity form and the gaussian critical temperature distribution 
used for the convolution. Panels C and D show the residual 
heat capacity after the convolved BCS heat capacity has been 
subtracted along with the convolved heat capacity, with Itac 
— 1.762 for panel C, and 5sc = 0.88 for panel D. The dashed 
lines are linear fits to the normal fermi liquid remaining above 
T 



within the cell. Above the bulk transition, Tco, a linear 
temperature dependence is seen, as expected for a normal 
Fermi liquid. We also note that in addition to the linear 
term in this region there appears to be a small additional 
constant contribution to the heat capacity (« 10/iJ/K). 
Golov and Pobell suggest that this contribution may arise 
from ordering in the amorphous solid ^He layer [ p^ . 

Since our main interest is in the heat capacity associ- 
ated with the '^He-aerogel superfluid transition, we shall 
subtract the bulk ^Hc heat capacity contribution. For 
this purpose we use the data of Greywall [ p6| and multi- 
ply by a suitable factor so as to produce a smooth con- 
tinuation of the remaining heat capacity data across the 
temperature of the bulk transition. The magnitude of the 
subtracted heat capacity indicates that a fraction (0.295) 
of the sample has the heat capacity of bulk ^He. This 
fraction is more than we can reasonably account for in 
terms of cracks and small volumes outside of the aerogel 
itself. Therefore we conclude that our sample of aerogel 
must contain a number of macroscopic pores. 

The heat capacity data remaining after the subtrac- 
tion of the bulk are shown in Fig. 3B. The heat capac- 
ity anomaly associated with the aerogel transition shows 
rounding but is clearly separated from the bulk transi- 
tion. The linear Fermi liquid region extends smoothly 
from temperatures above the bulk transition right to the 
transition in aerogel. In the region between the bulk 
transition and the aerogel Tc we find no indications of the 
anomalous behavior suggested by recent Grenoble NMR 
experiments 
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The ^He in aerogel effective mass, m*, can be deter- 
mined from the slope of the heat capacity in the Fermi- 
liquid region |16| |^0| after the fraction of bulk ■^He has 
been subtracted. From the slope of the data of Fig. 
SB above Tea, we find m*/m = 6.5, approximately 30% 
larger than the ratio for bulk superfluid '^He at this pres- 
sure. The apparent enhancement of m* cannot be as- 
cribed to a miscounting of the number of ^He atoms in 
the cell, and we believe that it must reflect a change in 
the excitations induced in the normal state by the pres- 
ence of aerogel. 

The simplest theoretical approach to the problem of 
superfluid ■^He in aerogel is that of the Homogeneous 
Scattering Model (HSM) |2|] , which treats the aerogel as 
homogeneous collection of scattering centers. The HSM 
predicts a reduction of the transition temperature com- 
parable to that seen in our experiments, as well as sup- 
pression of the order parameter. 

In order to assess the fraction of fluid contributing to 
the aerogel superfluid transition, we have fit the data to 
an interpolation of the BCS specific heat (Eq. 5) We 
take the weak coupling value of AC/C — 1.42, a transi- 
tion temperature. Tea, of 1.65 mK, t — T/Tca, the nor- 
mal fermi liquid heat capacity contribution, CN{Tea+), 



and the weak coupling value of S^c = 



_ A(0) 



1.764. 
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Vl(l + ^)e-"(^-)j, (5) 

The BCS specific heat is shown in figure 3B. Since the 
peak in the measured heat capacity is somewhat rounded, 
we fit it with a convolution of the BCS form with a gaus- 
sian distribution of critical temperatures, and display the 
residual heat capacity in figure 3C, along with the con- 
volved BCS heat capacity. The best fit is from a gaussian 
of width , and the gaussian is displayed in figure 

3B. The dashed line in figure 3C is a linear fit to the 
normal fermi liquid remaining above Tea- The residual 
heat capacity below Tea rises above this line, suggesting 
other contributions to the low temperature heat capac- 
ity. This might be expected as a consequence of the mag- 
netic ordering of the solid ^He substrate atoms, but this 
contribution is expected to be on the order of the tem- 
perature independent constant at our lowest temperature 
psf . Such an ordering of the solid layer is a feature ob- 
served in the earlier NMR aerogel experiments [p[ [^] ||] 
|§. As has been demonstrated by the NMR § [|[7 the 
influence of the magnetic solid '^He layer can be removed 
by the addition of a few layers of ^He on the substrate. 
This will be an interesting direction to pursue in future 
aerogel heat capacity experiments. 

The implied normal fermi liquid fraction is unchanged 
by choice of 5sc^ and remains at 0.30. If we apply the re- 
sults of Lawes and Parpia for the reduction of the energy 
gap for our T^a/Tco of 0.71 @, we find 5sc = 0.88. The 
residual heat capacity computed for this value is shown 
in figure 3D, and shows much less deviation from nor- 
mal fermi liquid heat capacity (drawn as a dotted line) 
than that of weak coupling energy gap value. At lowest 
temperatures, the residual rises above the normal fermi 
liquid, still suggesting an extra heat capacity contribu- 
tion, though much less than that implied from figure 3C. 

Somewhat surprisingly, we find that the fraction of the 
fiuid that appears to participate in the aerogel superfluid 
state is only a fraction, 0.30, of the total. Expressed in 
terms of the order parameter, one would conclude that 
there has been a reduction by a factor of over two under 
the conditions of our experiment. This result is in keep- 
ing with the order parameter reductions observed in the 
Cornell torsional oscillator experiments |^ |2^, and in 
agreement with the gapless superfluid model of Sharma 
and Sauls |2|. 

In conclusion, wc find that the behavior of "^He con- 
fined to aerogel shows normal Fermi liquid characteris- 
tics above the superfluid transitions. The enhancement 
of the effective mass along with the low temperature exci- 
tations contribution in this region can be further studied 
as a function of pressure in future research. While a size- 
able fraction of fluid in the cell exhibits bulk-like charac- 



teristics, a comparable fraction displays a BCS-like heat 
capacity centered at the reduced Tc, in conformity with 
expectations of a suppressed phase transition. Thus there 
appears to be a significant contribution (30%) of "bulk 
like" gapped excitations, and evidence for a gapless con- 
tribution to the specific heat in conformity to observa- 
tions in torsional pendulum measurements |0] . 
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